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Hydration of oxygen-deficient metal oxides causes filling of oxygen vacancies and formation 
of hydroxyl groups with interstitial structural protons, rotating around the oxygen in localized 
motion. Thermal activation from 500 to 800 K triggers delocalization of the protons by 
jumping to adjacent oxygen ions, constituting proton conductivity. We report quantitative 
analyses of proton and lattice dynamics by neutron-scattering data, which reveal the 
interaction of protons with the crystal lattice and proton-phonon coupling. The motion for the 
proton trapped in the elastic crystal field yields Eigen frequencies and coupling constants, 
which satisfy Holstein's polaron model for electrons and thus constitutes first experimental 
evidence for a proton polaron at high temperature. Proton jump rates follow a polaron model 
for cerium-oxygen and hydroxyl stretching modes, which are thus vehicles for proton 
conductivity. This confirms that the polaron mechanism is not restricted to electrons, but a 
universal charge carrier transport process. 
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T he proton is a structural element in hydrocarbons, metal 
hydrides and water. It is an elusive charge carrier 1 and the 
energy carrier in the hydrogen economy. As a structural 
element it is in a localized state; whereas as a mobile ion and as an 
energy carrier it should have a delocalized nature. 

Numerous metal oxides have proton-conducting properties 
and many rare earth transition metal oxides with perovskite 
structure are good proton conductors. Some of them are 
prospective solid electrolytes at temperatures from 500 to 
800 K 2 . The temperature-dependent proton conductivity 
is paralleled by substantial changes of the chemical and 
molecular structure 3 . Experimental studies show that the proton 
conduction is an Arrhenius-type thermal activated process, 
phenomenologically reminiscent of a polaron-type electronic 
conductivity 4 , which is observed in many 3d transition metal 
oxides with perovskite structure 5 . 

Therefore, a proton can be pictured as a trapped charge in 
an elastic crystal field and thus at least conceptually satisfy 
the definition of a polaron, which has been proposed and 
mathematically exercised specifically by theorists Samgin 6 and in 
more general terms by Krasnoholovets et al. 7 . Observations on 
the elastic properties made by experimenters Slodczyk et al. 8 with 
vibration spectroscopy made them speculate over the potential 
polaronic nature of the proton conductivity: ‘The question of the 
vibrational signature of isolated proton (for example, the ionic 
proton, a proton sharing its interaction with more than two 
acceptors) and its dynamic nature (proton gas, polaron and so 
on) is open’ 8 . Yet, such a proton polaron was not identified so far 
in hard inorganic matter such as metal oxides. 

Here we determine the proton jump times of hydrated yttrium 
substituted barium cerate with neutron-scattering and electro- 
analytical methods. The proton jump times follow a temperature 
dependence that fits the mathematical model of a polaron. 
Moreover, comparison with our vibration spectroscopy data 
suggest that it is the cerium-oxygen and hydroxyl stretching 
modes that propel the protons as charge carriers. 


Results 

Pressure-dependence of proton diffusion. We have investigated 
the proton dynamics in hydrated BaCeo.sYo.2O3 _ $ (BCY20) with 
quasi-elastic neutron scattering (QENS) in a specially designed 
cell for neutron-scattering experiments at high pressure and high 
temperature 9 . Figure 1 shows the QENS spectra of pressurized 
(P = 0.58 GPa) and non-pressurized (P = 0 GPa) hydrated BCY20 
at 720 and 770 K, respectively, at neutron momentum transfer 
(scattering vector) Q = 0.6 A - 1 and Q = 0.8 A “ l . 

In comparison with the neutron-scattering curves measured 
from a specimen under no pressure, the QENS line broadening in 
the high-pressure measurement is clearly noticeable. We can 
resolve and separate the QENS line broadening from the elastic 
scattering and the instrument resolution, as exercised in Fig. 2 in 
a previous work 10 . 

It is worthwhile to note that especially at 720 K, the spectra at 
around ±0.1meV are broader for the non-pressurized samples 
in the cell. This indicates that the QENS contribution from the 
non-pressurized hydrated BCY20 due to proton diffusivity 
disappears upon pressurizing the sample. This is a direct 
experimental proof that forced volume decrease causes suppres¬ 
sion of the QENS signature of the protons. In other words, when 
the crystal lattice is squeezed, the protons are slowed down. This 
is in line with the report by Azad et al 11 who reported two 
distinguishable cubic phases in BaZr 0 9 Y 0 JO2.95: one /Tphase with 
higher proton conductivity and larger unit cell and one a-phase 
with lower proton conductivity and smaller unit cell. At 770 K, 
the QENS difference between the spectra with and without 
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Figure 1 | Pressure effect on neutron spectra. QENS curves of hydrated 
BaCe 0 .8Yo.20 3 -<5 (BCY20) obtained in a high-pressure INCONEL cell at 
OGPa (red open data points) and at 0.58 GPa (blue filled data points) 
pressure at scattering vectors of Q = 0.6 A -1 at (a) 720 K and (b) 770 K, 
and 0.8 A _ 1 at (c) 720 K and (d) 770 K. The drawn lines are least square fit 
Gaussians accounting for the elastic contributions. 


pressure is smaller, because at sufficiently high temperature, 
770 K rather than 720 K, the protons resist the pressure and 
maintain higher diffusivity. 

For the separation of the QENS contribution from the elastic 
neutron-scattering curve, the structure factor S(Q,co) is subject to 
least-square fitting with the convolution of a Lorentzian 
quasielastic scattering function and the instrument resolution, 
as exercised in ref. 10. The Lorentzian full-width-half-maximum 
A(Q) is presented versus Q in Fig. 2 and constitutes a measure 
for the proton diffusivity. Smaller quasielastic line widths 
are distinguishable at 720 K for hydrated BCY20 measured at 
high pressure, which are around a factor of two smaller than 
the FHWM measured without pressure. 

Quantitative analysis of the width A(Q) from QENS spectra 
with the jump diffusion model by Chudley and Elliott 12 




(1) 


provides numerical results for the proton jump length / and time 
t between two jumps, which yields the proton diffusivity D via the 
relation l 2 = 6Dt. 


Activation energy for proton conduction. Figure 3 (top) 
shows the thus obtained diffusion coefficients D in Arrhenius 
representation vs reciprocal temperature. The slope drawn from the 
two diffusion coefficients obtained from the QENS of hydrated 
BCY20 powder in the INCONEL cell at no pressure yields an 
activation energy £ a of 0.35 eV (ref. 9). The corresponding data 
from the same powder under 0.58 GPa pressure yield a 
substantially higher value of E a =1.20eV. For comparison, the 
diffusivity coefficients from an independent QENS measurement 
on hydrated BCY20 sintered slabs were measured in a Pt 
container—the standard containment for QENS studies. The sin¬ 
tered slabs have obviously the lowest activation energy and an 
overall higher diffusion coefficient. As the diffusivity D can be 
expressed by the conductivity a via the Nernst-Einstein relation, 
we can transform the conductivities from electrochemical 
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Figure 2 | Accordance with Chudley-Elliott proton jump model. 

(a) P = 0 GPa. (b) P= 0.58 GPa. Q-dependence of the Lorentzian full-width- 
half-maximum A(Q) fit with the Chudley-Elliot model for hydrated 
BaCeo. 8 Yo. 2 O 3_,5 in the high-pressure INCONEL cell with pressure 
(closed data points) and without pressure (open data points), at 720 K 
(blue data points) and 770 K (red data points). The solid lines are the least 
square fits of the Chudley-Elliott model (equation ( 1 )) to the data points. 


impedance spectroscopy (EIS) and compare them with the QENS 
data: 


D = C 7- 


hT 

e 2 c H 


( 2 ) 


with Boltzmann constant /c B , temperature T, elementary charge e 
and concentration of protons c H per unit volume. In Fig. 3 (bot¬ 
tom), the diffusivities extracted from high-pressure EIS data from 
ref. 13 are presented for comparison with the QENS data. 

The change of the slope upon pressurizing the specimen in 
Fig. 3 (bottom) confirms also from the EIS data that the proton 
diffusion activation energy £ a increases with pressure. With 
increasing pressure on the specimen, the protons must overcome 
a higher activation barrier. This, together with our experimental 
findings for the proton conductivity at high pressure 9,13,14 , proves 
that the activation energy depends on the lattice strain parameter, 
which in turn is modulated by the elastic properties 15 . 

Higher activation energy in a smaller crystal lattice suggests 
that one single classical hopping process such as the Grotthuss 
mechanism 2 alone cannot explain our observation. Rather, 
cooperative processes such as lattice vibrations, that is, 
phonons, would assist in the proton transport. 

The mobility of protons has been rationalized as a consequence 
of phonon-assisted jumps 16,17 . Proton conductivity has been 
interpreted as a polaron-activated process 7 . Tomoyose et a/. 18 
found that the wave function of a proton between two adjacent 
oxygen ions in an AB0 3 -type perovskite would spread over the 
double well potential if the 0-0 separation was lower than 2.9 A, 
with the result that the proton would transfer between the energy 
minimum sites 18 . They derived the activation energy for this 
process from a small polaron model. 

In view of such striking analogy with the polaron-driven 
electric conductivity in some metal oxides, which manifests in the 
exponential T dependency of a (so-called Holstein mobility 
contribution 19 ), we considered a polaron-type charge transport 
also for our proton conductor. We follow in so far the spirit of 
Holstein 4,20,21 , and Yamashita and Kurosawa 22 whose work on 
polaron theory was inspired by the same observation 19 . 

The Debye temperature for BaCeY-oxide is 6 = 400 K (ref. 23). 
Holstein 4 shows that for T >1/2 6 the motion of the charge 
carrier is predominantly a diffusion process based on random 
jumps between neighbouring sites associated with Holstein’s 
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Figure 3 | Diffusion coefficients, (a) Diffusion coefficients D QENS of 
hydrated BCY20 obtained by QENS at OGPa (open red squares) and 
0.58 GPa (filled blue dots) in the INCONEL high-pressure cell, and at 
OGPa (open green diamonds) in a Pt container in the temperature range 
from 620 to 870 K. Data were obtained by least-square fits of A(Q) to the 
Chudley-Elliott model, (b) Comparison of diffusion coefficients obtained 
from QENS at 0 and 0.58 GPa (large symbols) and from EIS at 0.5-1.25 GPa 
(small symbols). 

off-diagonal transition. We are therefore considering the 
high-temperature region of Holstein’s small polaron theory. 


Proton polaron. The polaron concept was originally introduced 
by Landau for the description of charge transport in crystal 
lattices 24 , in particular for electrons, and developed further by 
other researchers 7,19,25 . Polaron theory was applied later to 
distinct different problems such as hydrogen diffusion in metal 26 , 
hydrogen-bonded systems such as ice 27,28 and perovskites 29-31 , 
for example. The proton with small ionic radius may have strong 
bonds with the lattice vibrations. When a proton is introduced 
into an ionic crystal and moves through the crystal lattice, 
because of their electrostatic interactions, the lattice ions respond 
to the proton motion by local lattice relaxation (polarization) 32 
as sketched in Fig. 4. The term ‘polaron’ implies that 
complete assessment of the problem would therefore include 
the effect of the ionic character on the pressure dependence of 
vibration frequency. 

Protons have an attractive interaction with negative charged 
O 2 - and a repulsive interaction with the positive charged 
Ce 4 + /Y 3 + in our BCY20 material. Rather than hopping within 
a static lattice, proton hopping can be understood as modulated 
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by the vibration of the Ce0 6 octahedra or by the O-H vibrations. 
Hence, the mobility of protons is a result of phonon-assisted 
jumps 6 . 


Phonon-modulated activation energy. With respect to the 
model depicted in Fig. 4, it is instructive to tentatively follow 
a classical approach by Wakamura 33,34 , who considers the 
vibration of a proton and an oxygen ion in the metal oxide 
host lattice as a damped harmonic oscillator. Accordingly, 
the relative displacement A u between a proton with mass ra H 
and an oxygen ion with mass m Q in an electric field with strength 

A M =l7 H -l7o= —(—-—) (3) 

ft>TO V m H7 H m oyo) 


Because of the inverse relation A u ~ l/co TO > higher optical 
phonon frequency co TO will mandate a smaller relative 
displacement. Assuming that proton hopping with activation 
energy £ a is assisted by a large vibration displacement of the 
oscillator, the activation energy scales parabolic with the optical 
phonon frequency: 


Fa OC 


1 

A u 2 


oc OJ 


TO 


(4) 


Therefore, we expect pressure dependence for the transversal 
optical (TO) and also longitudinal optical (LO) modes, and 
the different contribution of the vibration displacements for the 
high- and low-energy modes on the polaron transport. The 
former gives different pressure dependence on the LO modes and 
TO modes as we would expect from the magnitudes of the 
Griineisen parameters for the modes, and the latter gives different 
influence on the mobile proton. 

A change of the crystal unit cell volume will affect the phonon 
frequencies, the relationship of which is given by the Griineisen 
parameter y. The electronic polarization, which is in principle the 
optical dielectric tensor, affects the frequencies co TO and co LO of 
the TO modes and LO modes, respectively. Correspondingly, this 
yields different magnitudes on the Griineisen parameters for the 
LO and TO modes, y LO and y T0 . Wakamura 33 has demonstrated 
how different pressure dependencies arise for co TO and co lo . In an 
ionic crystal, y LO assumes a lower value than y TQ . As the polaron 
works stronger by the LO mode than by the TO mode, the LO 
mode must contribute dominantly to the activation energy F a . 
It may clarify the difference of parabolic and linear relations of 
vibration frequency on F a . Strong contribution of electronic 
polarizability on F a is obvious from the relationship between 
optical dielectric constant and F a values as seen in references 33,36 . 

Our EIS studies at high temperature and high pressure 13,14 
confirm the influence of pressure on the activation energy E a . 
At the molecular level, the relation between activation energy 
and mechanical pressure can be rationalized as the consequence 


of the shift in the prominent Raman mode v 3 , which represents 
the Ce-O stretching 15 . In other words, the frequency of certain 
phonon modes can affect the amplitude of the proton motion. At 
pressure, the spring constant (here not being constant anymore) 
of the Ce-O stretching bonds increases, which can result in 
smaller vibration amplitude of the protons. Pressure may keep the 
protons in the Ce0 6 octahedral cage and raise therefore the 
energy barrier, that is, F a for the proton transfer 9,37 . 

O-Ce-O bending modes were also observed in the Raman 
spectra of BCY at lower frequencies. The variation in frequencies 
for these modes upon pressure, however, was not as significant as 
the Ce-O stretching modes. In contrast, the pressure-dependent 
variation for these bending modes was suppressed upon 
hydration 37 . Therefore, although the O-Ce-O bending modes 
may contribute to proton transport via the vibration 
displacement 18 , the effect of these modes is not considered due 
to its insignificance here. 

Supplementary Fig. 1 (top panel) shows the proton conductiv¬ 
ity activation energies for pressures up to P = 2GPa, which we 
obtained with EIS from BZY10 (BaZro. 9 Yo. 1 O 3 ) and BCY20, and 
with QENS from BCY20. For each sample set, the activation 
energies vary quite linearly with pressure. Based on the linearity 
of the Raman shifts versus pressure 37 , we can express the pressure 
by wavenumbers. The bottom panel in Supplementary Fig. 1 
shows how £ a varies linearly with the Raman frequency (shift) 
from roughly 0.4 eV to around 0.8 eV in the measured pressure 
range. The flat grey line in Supplementary Fig. 1 (top) is the 
activation energy derived from equation (13) respectively from 
equation (14) in ref. 38, after we have expressed the square of the 
distance, u 2 > by pressure p : 

£ a =2/cT W 2 tan/ ! ( 8 ^ r ) (5) 


Thus, Samgin’s model predicts a linear relation between F a and p, 
which is in line with our experimental data in Supplementary 
Fig. 1 (top). We find no parabolic relationship in the investigated 
range of pressure and Raman shift. However, a virtually linear 
relationship of v and F a as was found by Wakamura 33,34,35 for a 
wide range of compounds for the OH stretching frequency 18 and 
was found in our study for the Ce-O stretching mode 
(Supplementary Fig. la). 

We also find no parabolic relationship between F a and co to . 
Instead, our investigation of the Raman shift under high pressure 
on hydrated BCY20 (ref. 37) produces a linear relationship 
between activation energy and phonon mode in line with the 
derivation by (refs 4,39) 


Fa 


h • ojq 
2n 


(6) 


For the elucidation of the difference of parabolic and linear 
relations of the vibration frequency on £ a , the effects of the ionic 
character or electronic polarizability would need to be considered. 
Further work needs to be done to study the difference between 
the pressure dependence of TO mode and LO mode for a solid 
conclusion. 


Proton jump time. Holstein 4 pictured the action of a small 
polaron as a random walk of a charged particle from site to site. 
The electron diffusivity is the product of the square of lattice 
distance and the total non-diagonal transition probability 4 . In 
analogy to Holstein’s model, Samgin 30 derived from the hopping 
mobility of his proposed proton polaron that under the thermal 
condition Fco ~ kT, the time t between two subsequent jumps 40 is 


t = (Am) *exp 


2 u 2 tanh 



(^) 
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Figure 5 | Polaron modelling of proton jump times, (a) Time i between 
two proton jumps obtained from QENS on hydrated BCY20, together with 
three least square fits of the data points to equation (7). Solid black line 
represents the Ce-0 stretching mode of 355cm -1 . Grey dotted and solid 
lines are O-H stretching mode of 3,560cm -1 (ref. 15) (b) Proton jump 
times obtained from EIS at various pressures using equation (7). Solid red 
line is the fit using equation (4) and the Ce-0 stretching mode (355cm -1 ) 
for the filled red squares measured at ambient pressure. The solid red line is 
the fit using equation (7) and the Ce-0 stretching mode (355cm -1 ) for 
the proton jump time derived from EIS data at ambient pressure. The fit 
works well for T>500K and deviates from low temperature data. No 
proper fitting to the EIS data can be obtained using the O-H stretching 
mode (3,650 cm -1 ). 


u is a coupling constant between proton and lattice vibration. 
With the experimental proton jump time t that we have directly 
determined by the least square fits from QENS data to the 
Chudley-Elliott model in ref. 9, we can verify the validity of the 
mathematical model for Samgin’s proposed proton polaron 
in ref. 30. Figure 5a shows how the time t between two proton 
jumps decreases systemically from 30 ps at 620 K to around 
7ps at around 870 K. The drawn lines are least-square fits 
to equation (7). This is experimental evidence that Samgin’s 
treatment of the Holstein polaron theory for the proton polaron 
is applicable to the proton transport studied here, that is, 
an experimental verification of Samgin’s proton polaron 
hypothesis and a validation and extension of the Holstein 
polaron theory. 

The proton jump time t can also be extracted from the 
diffusion coefficients obtained from EIS at various pressures. In 
the proton jump diffusion model, the following relation holds: 
/ 2 = 6Dt, where l is the proton jump length. We estimated 
the proton jump length from the 0-0 distance, obtained from 
high-pressure X-ray diffraction in reference 15 . We extracted t and 
plot in analogy to the proton jump time determined by QENS, as 


Table 1 | Proton-phonon coupling constant u obtained from 
least square fit to the lattice-assisted proton transfer model 
after Samgin, respectively, for the Ce-0 stretching and 
O-H stretching mode 15 . 


Method 

Assist phonon 
wavenumber (cm -j ) 

Factor A 
(cm ps -1 ) 

Proton-phonon 
coupling 
constant u 

Activation 
energy E a 
(eV) 

EIS 

QENS 

QENS 

QENS 

355 (Ce-0 stretching) 
355 (Ce-0 stretching) 
3,560 (O-H stretching) 
3,560 (O-H stretching) 

1.49 x10 12 
0.041 

7.6 x10 3 
2.83 x10 s 

5.29 

3.87 

3.23 (>700 K) 
4.00 (<700 K) 

0.64 ±0.03 
0.326 ±0.001 
1.29 ±0.21 
1.98 ±0.35 


EIS, electrochemical impedance spectroscopy; QENS, quasi-elastic neutron scattering. 
The activation energy E a is calculated according to ref. 30. 


shown in Fig. 5b. The solid red line is the least square fit of 
equation (7) using the Ce-O stretching mode (355cm -1 ) for 
the proton jump time derived from EIS data at no pressure. 
The fit, that is, the model reproduces the experimental data points 
well for T>500K, but deviates from low temperature data. 
No proper fit of the EIS data can be obtained using the 
O-H stretching mode (3,650cm -1 ). 

Table 1 summarizes the optimized fit parameters and the 
obtained proton-phonon coupling’ constant u (strictly speaking 
the coupling constant for the interaction of the proton with the 
lattice) for the respective vibrational mode. We first discuss a 
crystal phonon, the Ce-O stretching mode at 355 cm -1 , shown 
in the solid green line. The Ce-O stretching mode is the phonon 
mode, which shows the most significant change in frequency 
when pressure is applied 15 , that is, the response of this mode to 
the surrounding environment change is the strongest—the force 
constant of the Ce-O stretching vibration becomes ‘stiffer’ in a 
compact lattice. Therefore, when a proton-polaron is brought 
into the lattice, its interaction of the Ce-O stretching mode could 
also be strong. In fact, the wavenumber shift of 2 cm - 1 in this 
mode 15 is also the evidence that the phonon vibration changes 
due to the incorporation of protons. Samgin’s model suggests that 
the Ce-O stretching mode plays a significant role in the 
interaction with protons. 

The Raman feature at 3,560 cm -1 is a local vibrating mode 
and the signature of the O-H stretching. This mode appears 
in Raman spectra when the proton is incorporated in BCY20 
(ref. 15). The O-H stretching mode obviously promotes proton 
hopping for T>700K, as illustrated in Fig. 5a. Below this 
temperature, the solid grey fit line deviates from our experimental 
jump times. A further fit for T<7 00K is shown as a dashed grey 
line in Fig. 5a, resulting in different values for u and for a large £ a . 
This crossover in temperature coincides with other properties for 
the proton conducting perovskites, such as the discontinuity in 
the thermal expansion coefficient (second-order phase transition) 
and the change of the slope for the increase in conductivity versus 
the temperature 10 . 

This discontinuity in thermal expansion and proton diffusion 
has been observed for the same or for similar materials 41-44 . 
Literally, the thermal expansion coefficient decreases in the 
temperature range where we observe the onset of considerable 
proton diffusivity and conductivity, that is, here in the range 
650-700 K. Most noteworthily, the jump time t is correctly 
predicted for the Ce-O stretching mode by the proton polaron 
model (equation (7)), whereas the O-H stretching is validated 
only for T>700K. 

At low temperatures, the proton conductivity will follow 
a process, which is of quantum-mechanical nature 4,20 . 
Spahr et al. 40 reported that at low temperatures of 10-250 K, 
proton transport in rutile Ti0 2 is promoted by a transverse 
oscillation mode, which is coupled by O-H and O-D stretching 
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modes. In the case of the perovskite BaCeY-oxides, protons 
exhibit a localized rotational motion at low temperatures below 
500 K (ref. 45). 

Tomoyose et a/. 18 suggest a vital role of the bending mode 
of O-Ce-O in their polaron model. The difference from the 
Ce-O stretching modes as employed in this work is caused by 
the different proton diffusion path. Although the O-Ce-O 
bending modes can reduce the 0-0 distance within the same 
Ce0 6 octahedron, the Ce-O B lg and B 3g stretching modes as 
discussed in this work cause a ‘breathing’ motion octahedron 46 , 
which can lead to shorter distance between next-nearest 
neighbouring oxygen sites such as 03-03, as shown in Supple¬ 
mentary Fig. 2. The proton jump length obtained from QENS 
in ref. 9 suggests the existence of ‘inter-octahedra transfer’ in 
BaCeY-oxides, which is in line with an investigation using 
combined first principles and kinetic Monte Carlo simulations. 
Furthermore, the proton transfer between 03-03 sites happens 
only within the a-c crystal plane, which is in line with our finding 
that the phonon hardening occurs within the a-c plane 15 . 

Considering the motion of protons in the context of lattice 
dynamics, proton diffusion is caused by the cooperative atomic 
displacements. The phonons move mobile species towards the 
energy barrier. The polaron must overcome the potential barriers, 
which are posed by the lattice and modulated by its phonons. 
The protons are considered as hopping at sufficiently high 
temperatures from site to site due to exchange of energy and 
momentum with phonons 48 . The steps of phonon-assisted 
hydrogen diffusion in oxides are detailed in Fig. 6. The local 
relaxation of the oxygen ions induced by the proton is shown in 
steps (a) to (b). For long-range transport, a rotation step of the 
proton around the oxygen is necessary, that is, from (b) to (c). 
The new position of the proton in (c) results in another relaxation 
of the oxygen lattice, as shown in (d). The proton overcomes 
the energy barrier with thermal energy and hops to another 
oxygen site in (e). 

Effective mass of proton-polaron. Samgin 6,30 considers the 
proton polaron in a ceramic proton conductor a small polaron, 
because the range of interaction is smaller than the crystal unit 
cell, that is, a Holstein polaron. The polaron coupling constant u 
for the relevant Ce-O stretching mode is obtained from the least 
square fit of Samgin’s model to the experimental data for the 
proton jump times t from our QENS data and yields 
u — 3.87 ± 0.11. 

For intermediate polaron coupling with coupling constants 
<4.5, the following relation for the ratio of the polaron mass 
m* to the band mass m of the charge carrier is considered a good 
approximation 49 : 

— « 1 + - + 0.0236u 2 + 0.0014u 3 (8) 

m 6 

The thus obtained ratio between proton polaron mass and proton 
mass for the intermediate coupling is mVm = 2.0796 ± 0.0449. 
Within the weak coupling limit approximation, which includes 
not the cubic term u 3 , the ratio is mVm= 1.9985 ± 0.0382. 
The polaron mass is thus virtually the double proton mass. This 
outcome may point to the question whether the protons could 
actually move in pairs and thus constitute a proton bipolaron. 
A proton may attract the surrounding negative oxygen ions, 
moving the neighbouring oxygen ions towards the proton. The 
negative charge density in the vicinity of the proton is therefore 
increased. This negative charge could attract another proton, 
causing two protons to pair up. Interestingly, such bipolaron 
suggestion is consistent with our observation of inter-octahedra 
transfer in ref. 15 when the protons move in pairs, the nearest 
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Figure 6 | Stages of proton diffusion. Schematic representation of the 
phonon-assisted hydrogen diffusion process in a long range chain of 
oxygen, (a) Oxygen lattice with no protons, (b) Local relaxation of oxygen 
ions when a proton-polaron is incorporated, (c) Proton rotates around the 
oxygen to a new position, (d) The new proton position adjusts the oxygen 
lattice relaxation, (e) Proton hops to the next oxygen site. 


oxygen is already occupied by another proton; therefore, the 
mobile proton proceeds to the next neighbouring oxygen. 

Discussion 

In summary, our experimental observation that the proton 
jump times t are accurately described by Samgin’s proton 
polaron model lends credibility to the suggestion 6,7,8 that the 
proton-phonon coupling in an AB0 3 perovskite structure such as 
BCY20 can be considered the origin for a genuine polaron. In this 
case, the thermally activated Ce-O stretching mode is the 
dominant driving force for the coupling of the low-energy 
phonon mode with the proton diffusion, and thus the physical 
origin of the proton conductivity in ceramic proton conductors. 
The coupling strength of this system is within the range of 
intermediate coupling and yields a polaron mass of virtually 
exactly m* = 2. Although this result looks fortuitous, the 
suggestion that this could be an indication for a process, which 
is based on protons that move in pairs rather than individually, 
deserves attention and should warrant further experimental and 
theoretical investigation. 

Methods 

Sample preparation. The BaCe 0 .8Yo. 2 0 3 _ $ (BCY20) powder was prepared by 
solid-state reaction. Stoichiometric amounts of BaC0 3 (Sigma-Aldrich, 99%), 
Ce0 2 (Aldrich, 99.9%) and Y 2 0 3 (Stanford Materials, 99.9%) were ball-milled 
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in acetone with Zr0 2 balls of 10 mm diameter in a planetary mill ( 200 r.p.m.) 
for 1 h and calcined at 1,000 °C for 12 h twice in succession. The resulting powder 
was ball-milled again in a planetary mill and then axially pressed at 50 MPa for 
0.5 min into monoliths and sintered at 1,400 °C for 24 h. The obtained BCY20 has 
an average grain size of about 800 nm. For protonation, sintered monoliths or 
powder was heated to 500 °C in humid N 2 flow for 24 h, resulting in proton 
concentration equivalent to [OH] “ =3.8mol%. 

QENS at high pressure. QENS on hydrated BCY20 powder were conducted at 
720 and 770 K, at the pressure of 0 and 0.58 GPa in a high pressure cell, respec¬ 
tively 9 . QENS experiments were conducted at the neutron wavelength 2 = 5.75 A at 
the time-of-flight spectrometer FOCUS, SINQ Neutron Spallation Source 50,51 . The 
detailed cell pressurization, pressure determination procedure were described in 
ref. 9. QENS data reduction was made with the software package DAVE 52 . 

High-pressure impedance spectroscopy. In situ high-pressure EIS were 
measured on sintered powder at pressures ranging from 0.5-1.25 GPa in a 
piston-cylinder apparatus 53 from ambient temperature up to 440 °C. The EIS 
spectra were simulated in Software Novocontrol WinFit by an equivalent circuit 
model composed of three parallel RC series of the bulk, grain boundary and 
electrode responses 13 . 

Data availability. The data that support the findings of this study are available 
from the corresponding authors upon request. 
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